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Abstract 

Using a holographic dual model of QCD, we compute the pion electromagnetic form factor 
F-wiQ'^) in the spacelike momentum transfer region, as well as pion couplings to vector mesons 
Qpin)^^. Spontaneous and explicit chiral symmetry breaking are intrinsic features of this particular 
holographic model. We consider variants with both "hard-wall" and "soft-wall" infrared cutoffs, 
and find that the Fj^^Q'^) data tend to lie closer to the hard- wall model predictions, although both 
are too shallow for large Q^. By allowing the parameters of the soft-wall model (originally fixed 
by observables such as mp) to vary, one finds fits that tend to agree better with Ft^{Q'^). We also 
compute the pion charge radius (r^) for a variety of parameter choices, and use the values of 
Qpin)^^ and TTip"'* to observe the saturation of -F7r(0) by p poles. 
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I. INTRODUCTION 



Quantum chromodynamics, now in its fourth decade, has been known since its inception 
to stubbornly resist direct analytical solutions in its strong coupling regime. The modern 
period has seen numerous approaches developed to tackle this problem, sometimes by pro- 
viding a simplified picture of strong interactions and its states (e.g., quark potential models, 
chiral Lagrangians, quenched lattice calculations), or by working in energy or mass regimes 
where a key parameter may be assumed small (e.g., asymptotic freedom regime calcula- 
tions, operator-product expansions, heavy quark effective theory), or by studying features 
of quantum field theories that incorporate distinctly nonperturbative behavior (e.g., solitons, 
instantons). 

Perhaps the most interesting techniques are ones that stimulate advances in studies of 
QCD-like theories by making them more complicated by introducing additional degrees of 
freedom. In this category one includes the l/N^. expansion and, more recently, the grav- 
ity/gauge correspondence known by the names of the anti-de Sitter /conformal field theory 
(AdS/CFT) correspondence or the holographic dual approach l|. This proposed duality 
between strongly-coupled Yang-Mills theories and weakly-coupled gravity is exceptionally 
appealing because it implies a fundamental connection between gauge and string theories. 

The original example of this duality is given by A/" = 4 super- Yang Mills theory, which is 
conformal and therefore lacks asymptotic S'-matrix particle states. While this hardly seems 
like an auspicious starting point for modeling QCD-like theories and their rich hadronic 
spectra, a number of the original dual theories nonetheless possess such useful QCD-like 
properties as chiral symmetry breaking and confinement. Moreover, approximate conformal 
symmetry is a well-known property of QCD in the deep ultraviolet (UV) limit. 

In the holographic approach one begins with the 5-dimensional AdS metric, 

ds^ = gj^,[j^dx^^dx^ = —{rj^^dx'^dx'' — dz^) , (1) 

where t]^^ = diag(+, —,—,—) is distinguished from the full nontrivial 5D metric (^^^^ obtained 
from Eq. ([1]). The z (Liouville or "bulk") coordinate corresponds to an inverse energy scale 
{Q^l/z), in that the UV limit of QCD is represented by fields living on the AdS boundary 
z = (or, allowing for a UV cutoff, a small finite value z = e, a location called the "UV 
brane"). The gauge/gravity correspondence then states that every CFT operator 0{x) is 
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associated with a bulk field \l/(a:, z) uniquely determined by its value e) on the UV brane. 
The conformal symmetry is broken, introducing thereby a mass scale, by limiting the ability 
of the fields z) to penetrate deeply into the bulk, which corresponds to constraining the 
infrared (IR) behavior; this may be accomplished, for example, by imposing a hard cutoff 
and appropriate boundary conditions on z) at a value z = Zq (called the "IR brane") ^ 
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or by introducing a soft wall with an exponential decrease ~ e~'^ ^ in the action for large 
z jsl. The dimensionful parameters z^^ or k consequently serve the role of Aqcd- The 
Kaluza-Klein modes of the field z) then represent hadronic states of the same quantum 
numbers, producing towers of hadrons analogous to those arising in the original hadronic 
flux-tube string theories of the 1970's, which in turn are in close kinship with the original 
Regge theories of hadronic excitations. 

Fully exploiting the gauge/gravity correspondence to produce a model for real strong in- 
teraction physics — a method called "holographic QCD" or "AdS/QCD" — may be attempted 
either through a top-down approach starting with a particular string theory and choosing a 
background that (as mentioned above) naturally produces QCD-like properties, or a bottom- 
up approach starting with real QCD properties and using them to obtain constraints on 
viable dual gravity theories. In this paper we adopt the latter viewpoint. Work along these 
lines has become very popular in the past couple of years; especially well represented and 



close to this work in spirit are studies of hadronic spectra 



of hadrons in the presence of chiral symmetry breaking 
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Il7|,ll8|,ll9(. we are specifically interested in employing the 

formalism introduced in Ref. jd], in which chiral symmetry breaking is included directly in 
the Lagrangian, to obtain specific information on the form factor of the pion in this variant 
of AdS/QCD. 

According to the holographic correspondence, the global QCD symmetry of isospin asso- 
ciated with the two light quark fiavors is promoted to a gauged SU(2) symmetry respected 
by the bulk fields, and the problem reduces to one of constructing a 5D action containing all 
the fields of interest with the appropriate quantum numbers and Lorentz structures. Vary- 
ing the action with respect to the fields of interest generates their wave equations, which 
are solved as eigenvalue equations subject to appropriate boundary conditions to obtain the 
modes (mesons with particular masses and ^-dependent wave functions). Meson decay con- 
stants appear as values of the modes at the UV boundary (associated, as usual, with wave 
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functions at small distance scales) and form factors appear as convolution integrals over z 
of wave functions and currents. 



vector mesons (the p and its partners) 
4 1 using both hard- wall 1^ and soft- 



The behavior of wave functions and form factors o 

has been examined using the constructions of Ref. 

li 

wall [IS] IR boundary conditions. The former has the advantage of simplicity but produces 
the unphysical Regge trajectory ~ n^, while the latter produces the more phenomeno- 



logically realistic 
analysis of Refs. 
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behavior ~ [s]. A number of interesting results follow from the 



16| . including a distinctive pattern of vector meson dominance by 



the lowest-mass states and predictions of the p charge radius. However, the calculations 



of Refs. 
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IGj require only one dimensionful parameter, zq or k, respectively, and their 



eigenvalue equations for the vector mesons admit closed-form analytic solutions in terms of 
known functions (Bessel functions and Laguerre polynomials) . The pion sector as described 
in Ref. j3] also requires dimensionful parameters associated with both spontaneous (a) and 
explicit ivriq) chiral symmetry breaking, and the resultant wave equations cannot be solved 
in closed form for arbitrary values of the three dimensionful parameters. It is the purpose 
of this paper to examine both the analytic limiting cases of the pion interpolating field and 
numerical solutions for pion couplings and form factor, in both the hard- and soft-wall cases. 

In our numerical simulations we find that the data for the pion form factor Ft^{Q'^) 
lies closer to the hard-wall than the soft-wall model results. The three parameters of the 
hard-wall model were originally fit to mp, m^, and /,r, and we repeat the exercise for the 
parameters of the corresponding soft-wall model. One may adjust the parameters of either 
model to obtain a better fit to Ft^{Q'^), but at the cost of a poorer fit to at least one of these 
three observables. The F^^^Q^) data suggests an optimal model that incorporates features 
of both, but is closer to the hard-wall model. 

20 ] using the same formalism of Ref. j4| considers Ft^{Q'^) in the chiral 



A calculation 



limit {i.e., sets = 0) and focuses primarily on analytical behavior. We agree with their 
finding that the original hard- wall model overshoots Fj^lQ"^) data for large (but not with 
their choice of normalization for parameters such as /t^). 

This paper is organized as follows: In Sec. [Ill we recount and extend the formalism of 
Ref. relevant to our calculations and exhibit the analytically soluble limits to the equations 
of motion for the field containing the pion modes. Section IIIII gives expressions for the pion 
form factor and couplings in terms of AdS/QCD mode wave functions. Section HVl presents 
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the results of a number of numerical simulations of the pion electromagnetic form factor and 
couplings to vector mesons, and Sec. |V] summarizes our results and concludes. 



II. FORMALISM 



The full 5-dimensional action jj] used in this work reads 



a 



S = j dhe-''^ Tr + 3|X|2 - + F^jj 



(2) 



where g = \ det (^A/Af I obtained from the metric in Eq. ([T]), and e~* represents a background 
dilaton coupling, with $(-2) =0 in the hard- wall case and k'^z'^ in the soft-wall case. The 
Lagrangian within the braces of Eq. ([2]) is written in terms of a scalar field X and chiral gauge 
fields Al^r that enter through D^'X = d^'^X-iAf^^ X + iXAjl , Af^^j^ = A^^T with T being 
the generators of the gauged isospin symmetry, and F^^^ = d^^A^ j^— d^AfJ j^— i[AfJ j^, A^ pi\. 
The chiral gauge fields A'^\ are the holographic partners of the QCD operators ^ l^t^qj^ ^, 
while X"^ [more precisely (2/2;)X"^] is associated with q%qi, and therefore incorporates all 
chiral symmetry-breaking behavior. In the holographic dictionary of Ref. [4], the vacuum 
expectation value Xq of X (exact for the hard-wall model) is given by 

Xoiz) = ^Mz + ^j:z\ (3) 

where M = mgl and T, = al. represent explicit and spontaneous chiral symmetry breaking, 
respectively, and arise in the holographic recipe through the normalizable and nonnormaliz- 
able solutions for the bulk field X. Strictly speaking, Eq. ([3]) holds for the soft-wall model 
only for small values of 2; [3^; for large z the corresponding Xq should approach a constant, 
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but the e~'^ ^ background minimizes this distinction. We perform subsequent calculations 



using a background field that formally satisfies both limiting forms in Ref. 211]; however, 
as shown there, the best fits do not improve upon the naive soft-wall fits presented below. 
The pion field tt°' appearing through X = Xo exp{2iTT"-t°') is dimensionless and related to the 
canonically-normalized pion field vr" of chiral Lagrangians via '7r" = 7r°//,r, with /,r = 93 MeV. 
One now forms the polar and axial gauge fields V,A: V^^ = Aj^) and A*^ = 

- A^), in terms of which D^X = d^^X - i[V^^,X] - t{A^^,X}, F^^^ = d^'^V^ - 

QNyM_ ^ Q^A/ yTV] ^ [^Af ^7V])^ pMN ^ QM^N _ ^^^Af _ ^ QyAf ^ [^Af yJV]^^ 

S = yci5a;e-*V^Tr||Mp + 3|X|2-^(F^ + Fi)| . (4) 



This action represents the only terms quadratic in fields and containing two derivatives or 
less, and is therefore sufficient to obtain the free-field equations of motion for Vm, Am, 
and TT [generically z)]. It is convenient to work in an axial-like gauge, Vz{x,z) = 0, 
Az{x,z) = 0; the associated sources may then be expressed as divergences over only the 
usual four spacetime dimensions, d^V^ = (since isospin is conserved) and d'^A^^. A^ is 
further decomposed into a transverse (divergenceless) piece A^_[_ and a longitudinal piece 
ip\ = j_ + d^ip. The fields are chosen to satisfy Neumann conditions dz'^{x, z) = at 
z = zq in the hard-wall case and to give a vanishing contribution from the z ^ oo limit in 
the soft-wall case. 

Solving for the equations of motion of the fields z) Fourier transformed with respect 
to the 4D coordinates x yields 



1 + - — v 



0, 



z 



glv{zf^ 

H n (.TT — 



z-^ 



0, 



(5) 

(6) 
(7) 
(8) 



where v{z) = 2XQ{z)=mqZ + (yz'^ . Whereas Refs. [15|, ll6| obtain solutions of Eq. ([5]) (corre- 
sponding to the tower of p states), we focus instead on Eqs. ([7]) and ([8]), coupled equations 
corresponding to the vr and its excitations. The dimensionful chiral symmetry-breaking 
parameters niq and a render Eqs. ([7j) and ([8]) [as well as Eq. ([6])] more complicated than 
Eq. ([5]), which depends upon only one dimensionful parameter, zq or k. 

The gauge/gravity correspondence provides a method of determining the 5D gauge cou- 
pling g^. Applying the equation of motion Eq. ([5]) to the Fy portion of the action Eq. (jl]) 
leaves only the boundary term 



S=-7r^ d^'x V^dzV'^'' . (9) 

^fi'S J ^ z=e 

The significance of this quantity becomes clear when one resolves the vector field as 
V^{q,z) = V{q, z)V^{q), where V^{q) is the Fourier transform of the source of the vector 
current J^ = q'~i^t°'q at the UV boundary z = e, and V{q,z) (the "bulk-to-boundary propaga- 
tor") is normalized to V{q,e) = 1. Due to the isospin conservation constraint g^l^'^ = 0, one 



may replace V^V^"" with V^V^U^'^ 5'^'' and Yi'^" = 7]^^" —q^q" / ^ and then the usual quadratic 
variation of the action with respect to the source V produces the vector current two-point 
function: 



/ 



d^xe'^%j;{x)Jl{Q)) = 5^'U,,J:v{q'), (10) 
,2N _ e-^d,V{q,z) ^^^^ 



E 



from which one finds, matching to the QCD result for currents Jn normalized [16[ according 



to the prescription of {4] 



2 127r2 „ , . 

91 = ^-^4.^. (12) 



An analogous calculation in the axial sector relates the bulk-to-boundary propagator 
A{q, z) to the vr decay constant f.„: 

1 d,A{Q,z) 



f = — 

•In 2 



(13) 



The sets of normalizable eigenstates of Eqs. ©-([H]) form towers of hadrons of the corre- 
sponding quantum numbers. Since large is intrinsic to this procedure, the mesons have 
narrow widths and the spectral decompositions of self-energy functions such as are sums 
over poles: 

n=0 ^ " 

where Mn are the mass eigenvalues and /„ are the decay constants of vector modes ipniz) 
normalized according to 

dz— llJrn{z)lpn{z) = ^mn • (15) 

The coupled equations of motion Eqs. (IT])-® for the bulk-to-boundary propagators 
dz(p{q,z), dz^iqyz) can be combined to produce the decoupled and dimensionless Sturm- 
Liouville form 

[A(x) dM^)] + A(x) [q^ - /3ix)] y{x) = , (16) 

where, taking /i = 1/2:0 or k in the hard- and soft-wall cases, respectively, the dimensionless 
independent and dependent variables are x = jjiz and y{x) = [e'''^^^^^yx][d,j:(f{ixq,x/fi)], 
respectively. Furthermore, (f= q^ / and /3(x)= glv{fixY/x'^, or 

(3{x) = {rhq + axY , (17) 
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where mg= g^nig/ fi and a= g^cr/fi^ are dimensionless, and A(a;)= x/[P{x)e *(^/'')]. Equa- 
tion (IH]) then immediately gives a solution for the field dzir^q, z): 



da:TT(iJ,q, x/fj,) =q A{x) y{x) . 



(18) 



Equation ( 1161) is therefore a second-order ordinary differential equation in x with three 
dimensionless parameters, g^, rfig, and a. It appears not to admit a general closed-form 
solution in terms of well-known functions. However, one may, as in this paper, numerically 
solve the equation subject to physical constraints (fitting to m^, and m^). One may 
also explore limiting cases for various orderings of the parameters; since a enters Eq. ( IT6|l 
through the combination ax^, one may consider the behavior of ( JT6l) in the limits where 
various ratios of g^, (ax^)^, and are taken to be small. 

In the hard-wall case of $(x//i) = l, the analytically soluble cases are 
1) g^,m^ ^ (ax^)^ {z small or cr = 0) : 



y[x) 



AJq (^^q"^ - ml xj + BYq (^^ q^ - x j , or 



dzfiq, z) 



T - gi^K z + BYo [\ q^ - gim? z 



q^z 



1 9 2 

gim\z^Q L 



AJq \\J(f - glmlz) +BYo[Jq^- glm^ z 



(19) 



For extremely small or negative q^ (such that the arguments of the square roots become 
negative), the Bessel functions Jq and Yq of course analytically continue to modified Bessel 
functions. 

2) ^ (5"x^)^ ^ rhg {q^ large, z not extremely small) : 



y[x) 



x^ [AJ2{qx) + BY2{qx)] , or 



dMq.z) = -MUqz) + BY^iqz)] 



dzT^iq, z) 



^0 

^^.^^[AUqz) + BY2{qz)] . 



(20) 



The distinction between 1) and 2) arises from a noncommutativity of limits in A(x): If crx^ 
is taken small first, then A(x) -^x/rh,^ , while if m^^ is taken small first, then A(x) ^ 1/cr^x^. 
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(crx^)^ > > {z large) : 
y{x) = AAi' 



~ S\ 2/3' 



d,7r{q, z) 



3 \ 2/3 



~ 3\ 2/3' 



2 

+ BBi' 



or 



3 \ 2/3- 



gla'^z^zl 



AM 



(^) 



3\ 2/3- 



3\ 2/3- 



(21) 



Note that all the special functions appearing in this case are variants of Bessel functions, 
as seen repeatedly in previous papers that consider solutions to Eq. ([5]) for the hard-wall 
background. 

In the soft-wall case of $(x//i) = x^, the only analytically soluble cases turn out to 
have large g^, and the functions are variants on Kummer functions M(a, 6, z), t/(a, 6, z) (or 



22| : similar functions have 



equivalently, confluent hypergeometric or Whittaker functions) 

been seen for solutions to Eq. ([5]) for the soft-wall background [l6|. The analytically soluble 

cases are 

1) g^,m^ ^ {ax^)"^ {z small or cr = 0) : 



y{x) 



d,(fi{q,z) = K^ziAM 



e-^' {AM[l- _ rhl), 1, x^] +BU[l- i{f - mj), 1, x^] } , or 



1 — , 1, [f^z) 



d^q^z) = ^<AM 



A 2 ' 1' 



BU 
+ BU 



q^ - glml ' 
1 — , 1, (t^) 



4^2 



(22) 



2) q^ ^ {ax'^y ^ (g^ large, z not extremely small) : 



q 



y[x) = [AM(l-gV4, 3, x^) -gV4, 3, x^)] , or 

dMQ,z) = K^z^AM[l-q^/AK^, 3, (fi:;^)^] + B U[l - q^/ Ak^ , 3, (kz)'] } , 

9,7r(g,^) = ^^{AM[l-qy4K^3,{Kzy]+BU[l-q'/AK^3,{t^zf]} . (23) 

Numerous well-known recursion relations, distinct in form for M and U, may be used to 
reduce the arguments of the Kummer functions 22], but we opt to present expressions for 
which M and U have the same arguments. 
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III. FORM FACTOR EXPRESSIONS 



In the present calculations we are interested in the behavior of both form factors and the 3- 
point couplings Qmr-K [or 9p(^^)-k-k] between the nth vector state and the lowest eigenstate of the 
field TT. Of course, the n = case is the AdS/QCD version of Qp-m,. In order to identify these 
couplings, one must expand the action Eq. (jlj) to cubic order in fields. Since Eq. (jH]) relates 
the pion field to the longitudinal mode d^({) of A^, one must identify not only Vtiti terms, but 
also VAA and VAvr. Schematically, DX ~ dTx+T(dTx+0{-K^)+V-K+0{V-K^)+A+A'K+0{ATi^), 
Fv ~ dV + VV + AA, and Fa ~ dA + VA. X, the only field carrying SU(2) fundamental 
representation indices, must appear at least in pairs, while the dA term of F4 contains no 
longitudinal piece: d^{d'^ip) — d'^{d'^ip) = 0. The relevant terms then arise from the cross 
terms {dn){Vn) and iVTT){A) of \DX\^ and {dV){AA) of F^. One obtains the Vmr terms 

1 



^Ads = ^abc I d'^x I dz e" 



9b ^ 



(24) 



\2 

where the integration ranges over [0, Zq\ in the hard-wall case and [0, 00) in the soft-wall 
case 123! . Reference jj] uses this action (not including the dilaton coupling) to obtain the 
Vtiti couplings [Eq. (!33|) below], with the caveat that terms cubic in Fy^A have not been 
included. In fact, we now show that no such terms contribute to the Vtttt coupling. 

The field strength tensors Fy^ are antisymmetric and of opposite parities. Due to 
antisymmetry, only one Fy term (= Fy^ m Fy^ n Fy^ l) and one FyF\ term occurs, while 
the terms with an odd number of F^s are pseudoscalars. It is tempting to remove the parity 
distinction by forming the dual JF4 of F^-, but in 5D this object is a rank-3 tensor: 

^ JKL = q(^jklmn^^^^ ■ (25) 

Pairs of Levi-Civita tensors may always be converted into metric tensors using (5D versions 
of) the usual identities, so the only additional terms one might consider have a single JF4. 
But such terms cannot form scalars because they have an odd total number of Lorentz 
indices. Thus, only Fy and FyF^ need be considered. 

However, Fy terms that are linear in V contain at least 5 vr's, while the surviving terms 
of Fa contain at least one V, and therefore all terms in FyF]^ are at least quadratic in V. 
It follows that none of the F^ terms contribute to the Vtttt coupling. 
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Returning to the action Eq. (^^, a naive variation gives the 3-point correlator: 

mPi)J^'\q)J:i-p2)) = e'^'^Fiplplq') (pi + p^)" ^(2vr) - P2 + q) . (26) 

Again recalhng the narrowness of resonances, one may express the dynamical factor 
F{pI,P2, q^) in terms of transition form factors: 

17/ 2 2 2\ fnfkFnkijl ) /o'7\ 

where Fnk{(f') correspond to form factors ioi n-^k transitions. The pion form factor F^^lq"^) 
is then obtained as 

FAq') ^ Fuiq') = [ dz e'^ { ^[dM^)]' + Hz) - ^(z)]'] , (28) 

J Jtt I95Z Z ) 

an expression whose origin may be recognized in Eq. fl24l) . It is directly derived from that 
equation by factoring the 5D fields into products of the (dimensionless) 4D pion fields vr'*(g) 
and the bulk-to-boundary propagators tt{z) and ip{z). In order for the 7r'^(g) kinetic energy 
term to receive the standard canonical normalization |9^7r"(9^V", the integral in Eq. (IMl) 
[at = 0, where V{q,z) = 1] must equal unity, which fixes the normalization of Eq. (12^ . 
Note that the explicit factor in Eq. flM|) is part of the normalization and does not change 
the shape of F^(Q^). The pion is the ground-state solution to Eqs. ©-([H]) subject to the 
constraints at the large- 2; termini described in Sec. HIl The integral in Eq. fl28l) . setting 
V{q, z) = 1, is normalized to unity, giving a canonically-normalized kinetic energy term for 
the 4D pion field. Takin g a spacelike momentum transfer = —Q^ for the vector source V 



that solves Eq. ([5]) gives 



15| 



V{q,z) = J{Q,z)=Qz 



Ki{Qz) +Ii{Qz] 



(29) 



loiQzo) 

for the hard-wall case, while the corresponding expression for the soft-wall case is |16l . [1 

V{q, z) = J{Q, z) = r(l + Q^/Ak^) U[Q^IAk\ 0, [Kzf] . (30) 

Both of these solutions satisfy the boundary conditions V{q,e) = 1, V{0,z) = 1, as well as 
dzV{q, z) = for z = zq in the hard- wall case. In addition, for large z Eq. (150]) falls as 
^^2^-Q^/4^ Expressing instead V{q,z) for timelike momentum transfers gives 

°° f ip (z) 

Viq,z) = -g5j2 j^U2 - (^1) 



, g2 - M2 • 

n=l 
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Substituting this expression into Eq. (l28ll . one can represent the timehke pion form factor 
as a sum over vector meson poles: 

^^(^') = -Ei^' (32) 

n=l y " 

where gmrw is given by 

dzij^{z)e-''{^[dMz)f + '^Hz)-^{z)A . (33) 

Together, Eqs. fl28|) with Eqs. fl29|) or fl30|) and Eqs. fl32|) - fl33|) provide a complete expression 
for the pion form factor in all kinematic regions. 




IV. RESULTS 



In this section we present numerical predictions for QCD observables in both the hard- 



and soft-wall models. We perform the fit for the three hard-wall parameters Zq {z^ in 
mq, and a to the three observables rrip, m,r, and /,r- To the same observables we also fit 
the three parameters k, rriq, and a of the soft-wall model. In the hard-wall case the p wave 
functions [eigenfunctions of Eq. ([5])] are Bessel functions, with masses determined by zeroes 
of Jo{(1Zq)] hence, rrip = 70,1/^0 = 775.5 MeV, where 70,1 = 2.405 fixes 2:0 = 1/(322 MeV). One 
may then fit rriq and a to the experimental values of and /yr [which are constrained by 
the Gell-Mann-Oakes-Renner (GMOR) relation m^/^ = 2mq(T], yielding mq = 2.30 MeV and 
a= (326 MeV)'^ for the hard- wall model. In the soft- wall model the vector mass eigenvalues 
grow linearly in n [m?pn = 4(n + l)K^]; hence, the p (n = 0) mass fixes /t = mp/2 = 389 MeV. 
Note in particular that nip is fixed entirely by the value of /i = l/^o or k; iriq and a can then be 
fit to the experimental values of m.,^ and /tt, yielding for the soft-wall model mg = 1.45 MeV 
and a =(368 MeVf. 

Predictions for other QCD observables in both models are collected in Table HI Already 
one sees that the soft-wall model predicts fp to be much smaller than the experimental data, 
while the hard- wall model predicts a much closer value. Indeed, as noted in Ref. 3], the 
soft-wall model predicts the ratio m^/ fp to be exactly 2\/2t{ = 8.89, which differs dramatically 
from the experimental value 5.02±0.04 obtained from Table [H in comparison, the hard-wall 
model predicts this ratio to be v^7r7o,iJi (70,1) = 5.55. 

In Figs. [1] and [2] we plot the electromagnetic pion form factor Ft^{Q'^) in the spacelike 
region as obtained from both models. While results for both models are too shallow for all 
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Q"(GeV^) 



FIG. 1: Spacelike scaling behavior of Ft^{Q'^) as a function of = —q^. The continuous line is 
the prediction of the original hard-wall model. The dashed line is the prediction of the original 
soft-wall model with K = mp/2. The dash-dot line is the hard- wall model with (t= (254 MeV)^, and 
the dash-double-dot line is the soft-wall model with (t = (262 MeV)'^. The crosses are from a data 



compilation from CERN [2g] , the circles are from D 
the triangle is data from DESY 3Cll|, and the boxes 



SY, reanalyzed by Tadevosyan et al. 
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3l|], 



3ll | and diamonds [32^ are from Jefferson Lab. 



Older data in the range 3-10 GeV [33|] exist but have large uncertainties and are not plotted here. 

Q^, the hard- wall (solid line) prediction tends to lie closer to the experimental data than the 
soft-wall (dashed line) prediction. The hard- wall fit to Ft^{Q'^) can be improved for smaller 

by reducing 1/2:0 to about 255 MeV, at the expense of the vector meson parameters: 

1/2 

mp = 613 MeV and fp =260 MeV. In the soft-wall model, if one changes the value /t from 

388 MeV to 516 MeV in order to fit fp rather than mp (which becomes 1032 MeV), the 

result for F^(Q^) lies even further from the experimental data. On the other hand, if one 

lowers the value for k, which results in worse predictions for both fp and mp, then the fit to 

FniQ"^) lies closer to the experimental data; for example, k = 300 MeV (not plotted) gives 

the physical and values by taking a = (364 MeV)^ and nig = 1.70 MeV, but then 

mp = 600 MeV and /p^^ = 201 MeV. Lowering k to 255 MeV turns out to match F^iQ"^) data 

1/2 

somewhat better, but then mp = 510 MeV and fp =171 MeV. 
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Q^GeV^) 



FIG. 2: Spacelike scaling behavior of Q'^Ft^{Q'^) as a function of = —q^- The symbols are the 
same as in Fig. [TJ 

As noted above, once the vector parameters are determined by the value of /i = I/zq or 
K, the pion sector of the models determines best fit values for a and m^. One finds using 
Eq. (fT3l) that /tt depends mostly on a, and m^r is then fixed by choosing to satisfy the 
GMOR relation. Since, as is apparent in Figs. [T]and[2], the prediction from neither 

model is particularly good as increases, we consider the effect upon Ft^{Q'^) of varying a 
and TUq in both models. Empirically, both of the observables scale very close to the square 
root of the parameters (m^ oc rriq and oc a), but as one might expect, Ft,{Q'^) depends 
much more strongly upon a. Hence, if one allows a to float to fit the data for the 
precise fit to is spoiled. In particular, lowering the value for a to (254 MeV)^ in the 
hard-wall model but leaving 1/ zq = 322 MeV gives a much better fit to F^^IQ"^) (dash-dot 
line in Figs. [Hand [2]), but at the price of lowering the prediction of /,r to 64.2 MeV. In the 
soft-wall model, lowering the value for a to (262 MeV)^ and leaving k = 389 MeV gives a 
much better fit to F^^^Q"^) (dash- double- dot line in Figs. [Hand [2]), but at the price of lowering 
the prediction of /^r to 52.2 MeV. 

Using these form factor results for very low Q"^, one can extract the pion charge radius 
{rl) = -6dF^{Q^)/dQ\2=o. The experimental value (r^) = [0.672(8) fm]^ [2^ lies closer to 
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the original hard-wall [(r^) = (0.576 fm)^] than the soft-wall [(r^) = (0.494 fm)^] results, as a 
glance at Fig. [1] suggests. As remarked above, the hard- wall model fits the data better with 
a = (254 MeV)^ from which one finds (r^) = (0.645 fm)2, while setting cr= (262 MeV)^ for 
the soft-wall model gives (r^) = (0.600 fm)^. 

That the soft-wall model gives a shallower prediction for F^^^Q"^), and hence a smaller value 
for (r^) compared to the hard-wall prediction, is quite easy to explain numerically. First note 
that V{Q, z) in Eqs. fl29l) and fl30|) have quite similar z behaviors, except that V{Q, z) for the 
hard wall is cut off at z = zq. This fact alone allows for a greater contribution to the integral 
in Eq. ( l28l) in the soft-wall case. Moreover, the soft wall allows for more penetration of the 
Ti{z) and ^{z) fields into the bulk, as verified by our numerical simulations. So together 
they contribute more to the integration in Eq. fl25]) and give a higher value of F^((5^) at any 
particular value of than for the hard-wall model. 

Finally, we comment upon vector meson dominance and Fj^{q'^) in the timelike region. In 
the hard- wall case we find that -F7r(0) is essentially saturated by the first three p meson poles. 
Explicitly, from the values -Fn = 7o,n/[2oV^^I'^i(7o,n)|] and M„ = 7o,„/zo5 and using Eq. (l33l) 



to compute qit^t, = 2.3616, g2TTn = —0.8968, we find the first few contributions to ^,^(0) of 
0.8079-1-0.2830 — 0.0859 = 1.0050. However, the convergence is much slower in the soft-wall 
model. In the soft-wall case one has F„ = /t^(-\/2/vr)(n -|- 1)^/^ and M„ = 2K{n + lY^"^, and 
computes from Eq. fl33l) the couplings (^itttt = 3.3882, (727r7r = 2.9157, and (73^^ = 2.2946. With 
these values, we find the contributions to F^(0) of 0.3751+0.2696+0.1894+0.1291 = 0.9632. 
The next five terms are also positive and bring the sum of pole contributions to 1.138 before 



turning negative, a pattern that persists 
In contrast, results obtained in Refs. 



or several terms and brings -F7r(0) back close to 1. 
1^ for the form factor Fp{Q'^) show very different 



behavior, requiring fewer resonances for saturation in the soft-wall case. 

V. DISCUSSION AND CONCLUSIONS 

In this paper we have considered the problem of pion dynamical properties in holographic 
QCD, specifically in the context of models for which chiral symmetry breaking, both spon- 



taneous and explicit, is incorporated. In this way it differs from the recent work in Ref. |19| . 
which considers similar problems but contains no parameters analogous to a and rriq, and 
therefore uses a much simpler form for the pion bulk-to-boundary propagator. Nevertheless, 
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TABLE I: Hard- and soft-wall model predictions for QCD observables, the three model parameters 
in each case fit to m-,^, and nip (indicated by asterisks); all values except ^ptttt are in MeV. 



Observable 



Experiment 



Hard-wall 



Soft-wall 



m, 



ai 



.1/2 
Jp 

f 1/2 
Jai 



9pTl 



139.6±0.0004 [24] 
775.5±0.4 [24] 
1230±40 [24] 
92.4±0.35 [24] 
346.2±1.4 [25] 
433±13 [26, 27] 
6.03±0.07 [24] 



139.6* 
775.3* 
1358 
92.1* 
329 
463 
4.48 



139.6* 
777.4* 
1601 
87.0* 
261 
558 
3.33 



the basic result that Ft^{Q'^) is steeper in the hard-wall than the soft-wall model is common 
to both calculations. 

We also present explicit expressions for the pion bulk-to-boundary propagator dzTT{q,z) 
in all regimes where closed-form analytic solutions are possible. While not directly used in 
our numerical analysis of Ft^IQ"^), these results are useful for processes involving the tower 
of TT^") pseudoscalars. 

In our numerical studies of F^^IQ"^) we find that the naive hard- wall model appears to be 
more satisfactory than the soft-wall model over all regions of Q^, although the both models 
can be quantitatively improved by tweaking their parameters (at the expense of fits to other 
observables). Even before the F^^^Q"^) data is included, one can argue in favor or against 
either model based upon certain features; the soft-wall model has appropriate hnear Regge 
trajectories but poor agreement for the m'^p/fp ratio, while the hard-wall model is simpler 
and has the opposite behavior. The inclusion of the data for Ft^{Q'^) suggests an improved 
model with a semi-hard wall, such as provided by a quasi- "Saxon- Woods" background: 



^0 - 1 



(34) 



which has a drop-off at z = zq but falls off as e~^^^^ for large z, thus capturing features of 
both models and allowing for better predictions. The challenge, as always, is to predict the 
most observables with the fewest parameters. 
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